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ABSTRACT 

We present a pre-survey study of using Pan-STARRS high sampHng rate video mode 
guide star images to search for TNOs. Guide stars are primarily used by Pan-STARRS 
to compensate for image motion and hence improve the point spread function. With 
suitable selection of the guide stars within the Pan-STARRS 7 deg^ field of view, the 
lightcurves of these guide stars can also be used to search for occultations by TNOs. 
The best target stars for this purpose are stars with high signal-to-noise ratio (SNR) and 
small angular size. In order to do this, we compiled a catalog using the SNR calculated 
from stars with my < 13 mag in the Tycho2 catalog then cross matched these stars 
with the 2MASS catalog and estimated their angular sizes from (V — K) color. Wc also 
outlined a new detection method based on matched filter that is optimized to search 
for diffraction patterns in the lightcurves due to occultation by sub-kilometer TNOs. 
A detection threshold is set to compromise between real detections and false positives. 
Depending on the theoretical size distribution model used, we expect to find up to 
a hundred events during the three-year life time of the Pan-STARRS-1 project. The 
high sampling (30 Hz) of the project facilitates detections of small objects (as small as 
400 m), which are numerous according to power law size distribution, and thus allows 
us to verify various models and further constrain our understanding of the structure in 
the outer reach of the Solar System. We have tested the detection algorithm and the 
pipeline on a set of engineering data (taken at 10 Hz in stead of 30 Hz). No events were 
found within the engineering data, which is consistent with the small size of the data 
set and the theoretical models. Meanwhile, with a total of ^ 22 star-hours video mode 
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data < 10°), we are able to set an upper limit of A^(> 0.5 km) ~ 2.47 x 10^° deg-^ 
at 95% confidence limit. 



Introduction 



Since the first Kuiper Belt object (not including Pluto) was discovered in 1992 ( Jewitt Sz Luu 
19931, more than 1000 trans-Neptunian objects (TNOs) have been founcj^ A picture is emerging 
of the principal characteristics of the TNOs that is rich in information but leaves many unanswered 
questions. For example, a number of dynamical groups have been identified: classical, resonant. 



scattered disk (Luu et al. 1997), and more recently the extended scattered disk (Gladman et al. 
2002). What is the number and size distribution of the smaller objects of these dynamical groups? 
Is there an extension of the Kuiper Belt beyond 50 AU comprising bodies too small to have been 
detected in direct surveys? Meanwhile, could a possible close stellar encounter in the early history 
of the Solar System ( [Allen et al.|2001[ [Trujillo et al.|2001[ [Bernstein et al.||2004[ [Fuentes Holman 



2008) be responsible for the mass deficit and the depletion of larger objects (D > 150 km) beyond 
50 AU? 

The size distribution of these small TNOs provides important clues to the dynamical evolution 
of the early Solar System. Kilometer-size TNOs with 4% albedo are expected to have Mr > 30, 
which is still way below the detection limit of the largest ground-based telescopes. Yet the vast 
majority of the TNO population is beyond the limit of direct observation -for example, the Keck 
pencil beam survey ( [Chiang k Brown|[1999"| ) has a detection Hmit of Mr ^ 27.5. 

A number of authors ( Bailey|[l976 Rogues Sz Moncuquet[|2000 ) have suggested the possibility 
of detecting TNOs by stellar occultation. An occultation manifests itself as the shadow created 
by a TNO occulting a background star sweeping across the Earth, causing flux reduction in the 
lightcurve. Hence, as opposed to direct observation, stellar occultation by TNOs provides a unique 
way of detecting kilometer to sub-kilometer size objects in the foreseeable future. Recently a few 



groups have attempted to search for occultations by TNOs (Roques et al. 2006; Bickerton et al 
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detection has been confirmed yet. 

Projects committed to search for TNOs by occultation either do not have adequate time 
resolution, observe too few stars with good SNR or do not have enough observing time. These 
projects do a blind search for occultations, and these events are rare due to the reasons mentioned 
above. Therefore, the key to a successful detection is to have high sampling (> 20 Hz), high signal 
to noise ratio of background stars (SNR > 80), small stellar angular sizes (< 0.1 mas) and, most 
importantly, many star-hours of observations (> a few hundred thousands). 



^See 



http: / /www.cfa.harvard.edu/iau/lists/TNOs.html 



for a list of these objects. 
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The Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) is a project 
consisting of 4 telescopes that can cover over 6000 deg'^ per night or scan the whole visible sky from 
Hawaii (Svr) in a week to a detecting limit ~ 24**^ magnitude. Pan-STARRS prototype telescope 
(PSl) will monitor up to 60 guide stars in a very high sampling rate video mode. We have compiled 
a list which allows us to select guide stars from anywhere in the field; guide stars could be in any 
of the 64x64 OTCCDs (see ^ and read at ~ 30 Hz. Each field will have multiple choices of guide 
stars that were ranked based on the predicted event rates, which depend on their SNRs and angular 
sizes. By properly selecting guide stars, the Pan-STARRS video mode images would be ideal for 
searching TNOs. 

In the next section, we briefly describe the Pan-STARRS system. In ^ and ^ we discuss 
the diffraction profiles and present our detection algorithm. In § |5] we describe the prediction of 
stellar angular size using {V — K) color. In ^we estimate the number of expected events based on 
different number density estimations, sampling rates and stellar angular sizes. The methodology 
and compilation of the guide star list is in ^ In ^ we show and describe the quality of the 
engineering data obtained in fall 2008. Conclusions are in the final section. 



2. Pan-STARRS-1 system 

The full version of Pan-STARRS will contain four identical telescopes and camera systems. 



which is also known as Pan-STARRS-4 (Kaiser & Pan-STARRS Team 2002 1 . PSl is a single 



telescope and camera system built to test all requirements that are needed in the full version 
project. PSl telescope has an aperture of 1.8 m in diameter with 3-element wide field corrector 
that delivers a 7 deg^ flat field in the focal plane. Its 8 m focal length gives a plate scale of 38.5 
/im/arcsecond rendering a pixel resolution of 0.258 arcsecond/pixel. The telescope can monitor 
up to 240 guide stars at ~ 30 Hz (in reality, probably about 60 guide stars at 30Hz with high 
SNR) . High frequency image motion can be removed by on-chip fast guiding, where localized image 
motion can be calculated from a collection of nearby guide stars. Hence different image motions 
across the field of view can be corrected. This can only be achieved with the innovative design of an 
orthogonal transfer charge coupled device (OTCCD). The OTCCD allows on-chip charge shifting. 



in both rows and columns, which can compensate for the image motion in the focal plane (Tonry 



et al. 1997). The PSl Giga-Pixel Camera (GPCl) contains an 8x8 array, minus the 4 corners, of 
OTCCDs. For reasons of economy of manufacture, facilitating fast readout and other practical 
concerns, the OTCCDs each consist of an 8x8 Orthogonal Transfer Array (OTA) of independently- 
addressable cells of size 590x598 pixels. Hence GPCl comprises approximately 1.4 billion pixels. 
Typically, each OTA will have a single cell devoted to a guide star reading at high frequency. Of the 
590x598 pixels of this cell, only a small fraction around the guide star will be read, thus allowing 
a fast read out with relatively low read noise. 

PSl has four filters closely resembling SDSS filters (g,r,i,z), and a y-band filter. Their pass- 
bands, corresponding nii magnitudes and sky brightness are listed in Table [T] The major sources 
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Table 1:: Filter bandpass, mi (the AB magnitude that produces le /sec/pixel) and average sky 
brightness ii in magnitudes per arcsecond at Haleakala for PSl filters. 



Filter 


Bandpass (nm) 


mi mag 


fx mag/arcsec^ 


g 


405-550 


24.90 


21.90 


r 


552-689 


25.15 


20.85 


i 


691-815 


25.00 


20.15 


z 


815-915 


24.63 


19.26 


y 


967-1024 


23.03 


17.98 



of noise are: shot noise, background noise from the sky, read noise, radio frequency interference 
and dark current noise. The SNR for a given star can be estimated bjrl 



S/N = S/^ai + + ai^ + a^pj + 4 + ^7^ (1) 
where the various sources of noise are 

• Poisson noise, cxp = 0.5 x x t, where t is the exposure time and m is the apparent 
magnitude of the target star. 

• (sky) background noise, cig = x x t, where uj is the FWHM of the PSF in 
arcseconds, and /i is sky brightness in magnitude per square arcsecond; 

• read noise, = ^cj^ x Ax N^^^^, where A = 3.846 pixel/arcsec^ for 0.26 arcsecond pixels 
and A'rcad is the read noise of the detector. 

• noise from radio frequency interference, cr^pi ~ f "-^^ xAx A^p^pj, where N^pj is radio frequency 
interference noise. 

• noise from dark current, cr^ = jlo'^ x A x d x t, with d counts/sec 

Based on the characteristic noise mentioned above, we calculated the SNR versus magni- 
tude with different filters under 30 Hz sampling, 27e~ read-out noise, 1 arcsecond FWHM, d = 
5 counts/sec and no radio frequency interference. Figure [T] shows the relation between magnitude 
and SNR for all filters. Except for the y filter, stars with my < 11 mag can achieve SNR > 80 with 
30 Hz sampling. 



^We took the numbers from PSl document PSDC-230-002-00 (Mission Concept Definition http://pan- 
starrs.ifa.hawaii.edu/project/PSDC/PSDC-2xx/ I for the characteristic noise for all possible sources. 
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3. Diffraction Profile of Stellar Occultation 

The differential size distribution dN/dD oc implies that large objects are small in number. 
Therefore, the odds of serendipitous occultation by a large object is very small. Instead, we stand 
a better chance looking for occultations by objects in the small end of the size spectrum. There 
are two reasons for this: kilometer to sub-kilometer TNOs are far more numerous than big TNOs 
as the size distribution suggests, and the effective diameter of the occultation shadow is enlarged 
due to the diffraction effect, which effectively increases the event rates. 

In diffraction theory, there are two regimes delimited by the Fresnel number F = D'^/{rX), 
where D is the size of the TNO, r is the distance from the object to the observer, and A is the 
observing wavelength: the Fresnel diffraction (near field) where F and Fraunhofer diffraction 
(far field) where F \. As we increase the size of the TNO or decrease the distance to the 
object, we pass from the Fraunhofer regime to the Fresnel regime. The characteristic size, known 
as Fresnel scale F 2 km at 43 AU, 600 nm), is defined by \/rA and can be used to identify the 
two regimes. The regime of F ^ 1 is called the geometric regime, however for the scope of this 
work the probability of having such an occultation is insignificant. 

The diffraction profile of stellar occultation for a point source can be calculated by the Lommel 
function and the profile for stars with finite size can be obtained by convolving the point source 



profile with the stellar disk (Roques et al. 1987). The profile is a function of TNO distance, size 



and shape, stellar angular size and impact parameter (defined as the distance between the line of 



sight of the star and the object). Nihei et al. (20071 showed that the width of a diffraction profile 



H (the diameter of the first Airy ring) with zero impact parameter is defined as: 

2 

H{\r,D,e,) = 2y{^JJj2F)l +{D/2)^' +re,, (2) 

where 6^ is the angular size of occulted star. Figure[2] shows diffraction profiles for various TNO 
sizes at 43 AU with zero impact parameter at opposition. As the size of TNOs becomes smaller 
than the Fresnel scale, the effective diameter remains almost the same. In other words, for objects 
smaller than the characteristic length defined by Fresnel scale, the diffraction effect will make the 
effective diameter independent of the size of the object. 

For non-point sources, due to the fringe superposition, stars with large stellar angular size 
would smear out the diffraction profile and reduce the depth of occultation as shown in Figure|3j 
the diffraction profiles of 600 m size TNO at 43 AU for various stellar angular sizes. 

For more distant objects, the Fresnel scale gets larger, and thus the same size objects are in 
the Fraunhofer regime. Furthermore, the projected size of the background star gets larger, thus 
a distant occultation event would have longer event duration but shallower occultation depth. In 
Figure|4j as the distance to the objects increases, the effective diameter becomes larger, and the 
depth of the flux reduction gets shallower. 



The diffraction profile of spherical objects have very distinct Poisson peak which can not be 
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Fig. 1. — : SNR plot for PSl filters with 27e~ read out noise, 1 arcsec FWHM, d = 5 counts/sec 
and 30 Hz sampling. Sky brightness and mi values for each filter are taken from Table [T] Note 
that, except for the y filter, stars with my < 11 mag can achieve SNR > 80 with 30 Hz sampling. 
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Fig. 2. — : Diffraction profiles for various TNO sizes at 43 AU and zero impact parameter at 
opposition. As the size of TNOs becomes smaller than the Fresnel scale, the effective diameter 
remains almost the same. In other words, for objects smaller than the Fresnel scale, diffraction 
effect will make the effective diameter larger than the geometric diameter. The Fresnel scale at 
43 AU is about 2 km, therefore in this case, 300 m, 700 m and 1000 m objects all have the same 
effective diameter. 
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Fig. 3. — : Diffraction profiles of a 600 m size TNO at 43 AU and zero impact parameter at 
opposition for various stellar angular sizes. Due to the fringe superposition, a star with large stellar 
angular size would smear out the diffraction profile and reduce the occultation depth. As it can be 
seen from the differences between 0.1 and 0.01 mas stars, the diffraction profile is very distinct for 
0.01 mas star and becomes shallower when the stellar angular size gets larger. 
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Fig. 4. : Diffraction profile for 1 km size TNO at various distance with zero impact parameter at 
opposition. The effective diameter becomes larger and the depth of the fiux reduction gets shallower 
as the distance to the TNO increases. 
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observed in irregular shape objects. However, as the objects become smaller than the Fresnel scale, 



the central Poisson peak exists for all shapes of objects (Roques et al. 1987). 



Finally Figure [5] shows the diffraction pattern for various impact parameters. When the objects 
are smaller than the Fresnel scale, the Poisson peak is prominent for 6 = 0. The peak falls gradually, 
as b increases, to a trough and then the intensity reverses (brightening event) as it reaches the first 
Airy ring. 

One can notice that in the Fraunhofer regime, for example D = 500 m at 43 AU, the shape and 
size do not change the width of the diffraction profile, only the depth of the occultation is affected 
by the object size and impact parameters. We based our detection method on this property as we 
explain in the following section. 



4. Detection Method 



When the shadow of a stellar occultation sweeps across a telescope, high sampling rate systems 
like PSl video mode will translate the shadow with diffraction profile into several time-sampled 
lightcurve points. In order to look for signatures of occultation, a detection method should be able 
to probe these lightcurves for diffraction patterns. There are far more numerous TNOs smaller 
than the Fresnel scale, therefore we target our search at objects with D < F. In this regime, H 
does not change as a function of the object size, shape or stellar spectral type. On the other hand, 
the depth of the occultation depends on the object size and stellar angular size. Also, the shape of 
the diffraction profile depends on the impact parameter. 

Based on these conditions, we have designed our detection method to search for occultations 
by scaling pre-defined TNO diffraction profiles (i.e. templates) to match patterns in the lightcurves. 
For a given star, with known stellar angular size 6^ and SNR S, let ^(r; b, {6^, S}) be a normalized 
template with mean zero as a function of time r and impact parameter b. We then try to find the 
best matched template by maximizing the following measure, 

T=t+'W 

M{t)= [ag{T;b,{9.,S})-p{T)f , (3) 

T=t — W 

with respect to the scaling factor a, where p is the local nominal baseline subtracted lightcurve 
within a window [—w : +w]. The size of the window n can be found analytically as sH/v^ei, where 
s is the sampling rate, v^ei = Ve (cos (p — \Jt^^) is the relative velocity with Ve is the orbital speed 
of the Earth, cj) is the opposition angle and r^^ is the distance of a TNO to the observer in units 
of AU. Finding the scaling parameter a can be done analytically as: 

^(u tn c^^ _ Ep{r)9{T;b,{e^,S}) 

l^g{T;b,{9^,S}y 

To consider the effect of b, we simply maximize a over b, 
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a{{e,,S})=max[a{b,{e,,S})] . (5) 

b 

In practice, for a given 0^, and S, we built a set of templates with impact parameter values, 
b = {0.00H,0.01H,0.02H, . . . ,0.50H}, and scan all templates across the lightcurve to find the 
template that maximizes a{b, {9^, S}). In order to speed up the computation, we actually perform 
a two-pass procedure. The first pass involves scanning corresponding templates from b = O.OH to 
0.5H with Ab = O.IH to find the region that brackets the maximum a. Then within this region, 
we apply a second pass that uses templates with finer impact parameter step, Ab = 0.01/7, and 
again find the maximum a value within this region. The template that gives the largest a value is 
the best fit between the scaled template and a pattern in the lightcurve. At the end of this process, 
we have the same number of optimized a's as the number of points in all lightcurve^ In order to 
decide if any of these a values represent real events or false positives, we determined a threshold 
value at{0^,,S), calculated from the null hypothesis distribution. For this purpose, we simulated 
lightcurves with no events, then scanned templates through these lightcurves, optimized for a's, 
built the distribution, set the p-value based on the desired false positives and finally determined 
the at{6i,,S) accordingly. 

This distribution can also be estimated semi-analytically, thereby reducing the computational 
cost. Finding an optimized parameter using a least square metric while minimizing over another 



parameter is equivalent to a regression with one nuisance parameter (Rice 2007 1 . The cumulative 
distribution of such a system is known to follow the Gumbel distribution or type I extreme value 
distribution ( Gumbel]|1943 1, 



G{X <x) = e---^'^-'^' , (6) 

where v and ^ are fitted parameters (scale and location parameters). 

The left panel of Figure[6] compares the null hypothesis cumulative distribution of ci's, for 
SNR=100 and 9^ = 0.01 mas, built from simulations (solid points) with the fitted Gumbel distri- 
bution (solid line). As can be seen from the figure, the fitting for cumulative distribution < 0.8 
(d < 0.2) was good. But the fitted part of the tail {a > 0.2) of the distribution, which is the 
part we care about the most, did not fit well. Therefore, we separately fitted the tail part of the 
cumulative distribution, and the results are shown in the inset of left panel of the same figure. 
Note that the difference in the tail between the null hypothesis distribution and the fitted Gumbel 
distribution is now very small. The right panel of Figure[6] shows the number of false positives 
predicted by Gumbel distribution versus the number of false positives observed in the simulation. 
As can be seen, at high threshold (bottom-left corner) the number of predicted false positives is 
more than the number of observed false positives. This is due to the fact that the at{0-^, S) set by 
fitted Gumbel distribution is actually higher than that of the true null hypothesis distribution. As 



^There are actually fewer a's than number of points because at the beginning and end of the lightcurves there are 
not enough points to define the local nominal baseline. 
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a result, even though the threshold will be set a bit higher, at(^*,5') can be still calculated based 
on the Gumbel distribution to ensure one false positive in the lifetime of the project. In reality, we 
may get fewer than one false positive, as we demonstrated above. 

The next step is to search in the lightcurves for any a values that are larger than the at{6^, S). 
Note that the at{9i,, S) is 9^, and S dependent, thus the process of simulating the null hypothesis 
distribution and fitting for u and fi has to be repeated for each individual star. In order to 
alleviate this problem, we have constructed a grid of fitted at(^*) 'S')'s values for various SNRs and 
angular sizes a priory and tabulated these values into a database. For a given star, the threshold 
is approximated with at{9'_^, S') where 9'^ and S' are the closest pre-calculated values to 9^ and S. 

Note that in the presence of auto-correlations (trends) in the lightcurves the Gumbel distri- 
bution may not best describe the data. We have not yet investigated how well this approximation 
performs when trend exists in the lightcurves. Therefore, in the event detection pipeline, we apply 
a high-pass filter and de-trend the lightcurves to avoid this difficulty (see 

To estimate the detection efficiency as a function of TNO-size for PSl guide stars with various 
SNRs and angular sizes, we implanted objects with various sizes into simulated lightcurves and 
calculated the number of a values for each object size that are larger than the threshold value 
at{9i,,S). Figure [7] left panel shows the detection efficiency for different stellar angular sizes as a 
function of object size. As the plot suggests, PSl can detect objects as small as ~ 400 m. Figure[7] 
right panel shows the detection efficiency as a function of SNR and 9^, for 700 m obejcts; the 
detection efficiency is higher for higher SNRs and smaller stellar angular sizes. Therefore the guide 
stars selection (see Q should be based on this fact: stars with higher SNR and smaller stellar 
angular size are best candidates for searching stellar occultation in the lightcurves. 

To evaluate whether our detection process is computational feasible, we have estimated the to- 
tal float point operations (FLOP) needed to process all data for one night to be ~ 10^^ (this is a con- 
servative estimate). To save computational time, we firstly tabulate the values of ^ g{T; b, {0*, S*})^ 
in Equation|4| The two-pass procedure to search for matched template needs about 25 FLOP. Each 
template has 13 data points, therefore for each data point in the lightcurve, the computational cost 
is about 325 FLOPS and for each night, there will be a total of ~ 5 x 10^ data points assuming 
30IIz sampling and 60 guide stars. With a 3G Hz CPU, the whole detection process can be done 
in few minutes. 



Bickerton et al. (2008) has demonstrated a similar detection method using convolution of 
templates with lightcurves. The advantage of convolution done in Fourier space is that it is com- 
putationally faster. Nevertheless, our computational cost is not significant as we have pointed out 
above, and therefore we choose to do the template matching in real space in order to avoid problems 
of aliasing and for better control and understanding of the underlying distribution. 
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Fig. 5. — : Diffraction profile for 500 m size TNO at 43 AU with various impact parameters at 
opposition. The central Poisson peak drops as the impact parameter increases until it reaches the 
first Airy ring where the intensity goes up again. 
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Fig. 6. — : Left panel shows the null hypothesis distribution of a values (shown for SNR=100 
= 0.01 mas with dots) and fitted Gumbel distribution (solid line). From a separate fitting to 
the tail part (inset), the difference between the null hypothesis distribution from simulation and 
the fitted Gumbel distribution is very small. The right panel shows the number of false positives 
predicted by Gumbel distribution versus the number of false positives observed in the simulation. 
The diagonal line is to show the difference between predicted and observed values. As can be seen 
that at high threshold (bottom- left corner), the number of predicted false positives is more than 
the number of observed false positives. 



5. Estimate of The Stellar Angular Size 



As mentioned before, the stellar angular size is important in choosing the guide stars. For this 
reason we need to be able to estimate these angular sizes ahead of time in order to select guide 
stars that will yield the highest number of events. 

Various methods have been used to measure stellar angular sizes -for example, lunar occulta- 
tion, stellar interferometry, speckle interferometry, and other photometric and spectrophotometric 



methods (see Scholz (1997) and reference therein). However, these measurements from the litera- 



tures are limited to relatively few bright stars. An alternative way to select the most appropriate 
guide stars from all stars visible to PSl Svr sky, is to estimate the stellar angular sizes using their 



colors. For example. Van Belle] (1999) fitted a few stars with known zero magnitude angular sizes 
{6v=o and 6b=o) versus their {V — K) and {B — K) colors to set up the relation between the color 
indices and their angular sizes. Accordingly, we used the formulae in his work in combination with 
the Tycho20 catalog ( |H0g et al.||2000[ ) for F-band magnitudes, and the 2MASS0catalog ( [Skrutskie 
et al.|[2006| for i^T-band magnitudes. We first extracted all stars with my < 13 mag and mx < 16 
mag above the PSl southern declination limit 6 > —30° and cross matched stars from each catalog 
by looking for pairs within one arcsecond radius and calculated their corresponding (V — K) color. 
We then refitted recent measurements of angular sizes from the Catalog of Apparent Diameter and 



Absolute Radii of Stars (CADARS) compiled by Pasinetti Fracassini et al. (2001) with Van Belle 



(1999)'s relations. We obtained the following fitted relation: 



= ]^q[0.453±0. 003+0. 246±0.005(y-ft:)-0.2y] 



for main sequence stars and 



-|^q[0.407±0. 007+0. 238±0.009(y-ii:)-0.2y] 



(7) 



(8) 



for giants. The fitting results are shown in Figure[8] and as can be seen, the predicted values 
matched well with the measured ones for both main sequence and giant stars. 

Due to the fact that stellar angular sizes depend not only on the spectral type but also on the 
luminosity class, the next question is how to reckon the luminosity class for every matched star 
when there is no such information. To do so we relied on the {V — K) color as a way to (roughly) 
distinguish between main sequence and giant stars. On the left panel of Figure [9] we show that the 
distribution of iy — K) from all matched stars in PSl 37r sky has two peaks, suggesting two groups 
of stars - most likely main sequence and giants. We used the luminosity class information provided 
in CADARS to determine the boundaries between these two classes (shown in Figure |9] middle 
panel). Based on this, all stars with [V — K) < 1.85 were considered to be the main sequence 
stars and all stars with (V — K) > 2.0 to be giants. We abstained from classification for stars 



''http://archive.eso.org/ASTROM/TYC-2/readme.htm 
^ http : / / www. ipac . caltech.edu /2mass/ 
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Fig. 7. — : The left panel is the efficiency as a function of TNO size for SNR=100 and various stellar 
angular sizes. The right panel is the efficiency contour as function of SNR and stellar angular size 
for 700 m objects. Stars with higher SNR and smaller stellar angular size have higher detection 
efficiency and are best candidates for searching stellar occultation in the lightcurves. 




Fig. 8. — : Angular sizes predicted from Equation[7] and Equation[8| The predicted values matched 
well with both main sequence (left panel) and giant (right panel) stars from CADARS. 
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between these boundaries. Note there are stih a few giant stars with small (V — K) (see middle 
panel of Figure|9]). This is a selection bias in the catalog due to the fact that these are bright blue 
O or B giant stars with multiple measurements. In reality, O and B giant stars are rare and can 



be removed from the final catalog using known O and B star catalogs (Reed 2003; Maiz-Apellaniz 



et al. 2004). 



Using the {V — K) cut, Equation[7] and Equation[8j we predicted angular size for CADARS 
stars and estimated the rms error to be 21.7% for main sequence and 20.6% for giants. 

Finally we estimated the angular sizes for all matched stars visible to PSl Svr sky and the re- 
sulting distributions are shown in Figure[9] (right panel). Main sequence stars peak around 0.02 mas, 
and giants peak around 0.05 mas where the overall distribution (solid-line) peaks around 0.02 mas. 



6. Event Rate Based on Various Models 

In this section we estimate the number of expected events for the three year project lifetime 
based on different models. The comparison of the number of expected events with the actual 
observed events will allow for discrimination between the different models. 

One can estimate the number of expected events as a combination the total solid angle, detec- 
tion efficiency and differential size distribution. It follows that: 



Advent - 5Z ^("^ife ' ^ / 



°° dN 

X 



dD 



k jk 

H{X,r,D,ek) e{X,D,ek,Sk)dD , (9) 

where n is total number of guide stars, jk is the jth exposure for guide star k, v{(pj^,r) is the 
apparent motion of the TNO for guide star k and exposure jk, ^j,. is the exposure time, dN/dD 
is the differential size distribution, H{X, r, D, Ok) is the effective diameter, and e(A, D, 6k, Sk) is the 
detection efficiency for PSl. The apparent motion v{(pj^,r) and effective diameter Hk{X,r, D,6k) 
can be calculated with = 0° (PSl will observe mainly close to opposition) for all stars and all 
pointings, r = 43 AU and A's corresponding to the particular filters. The total number of exposures 
Ujf, for each star k will be reduced according to the time that PSl is pointing within the ecliptic. If 
we assume that TNOs are distributed uniformly within < 6° along ecliptic, the total time PSl 
will spend near ecliptic is then sin(/?) r ~ 0.10 r, where r is the total time of observation in seconds 
for a three-year lifetime. The detection efficiency for PSl has been investigated as a function of 9k 
and Sk in f|4] The last unknown in the event rate calculation is the model dependent differential 
size distribution. For this, we adapted three differential size distributions. First, we used a double 



power-law distribution described in Bernstein et al. (2004 1 (hereafter B04) with bright-end slope 



ai = 0.88 and faint-end slope a2 = 0.32. The conversion from R-magnitude to size of the TNO 
was done by using 4% albedo and cf) 0° ( Russell||1916 1. The cumulative size distribution is shown 



as B04 line in Figure 10 Second, we use the models by Pan &; Sari (2005) (hereafter P&:S05) 
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with three different break diameters at 20, 40 and 80 km and g = 3 also shown in FigurefTOj The 



third model is from the upper limit set by the TAOS project (Zhang et al. 2008) (hereafter Z08). 
Anchored at 28 km, TAOS collaboration reported an upper limit on the cumulative surface density 
for objects larger than 0.5 km shown as solid circle (anchored point) and solid square (upper limit) 



also in Figure 10 Other models, such as Kenyon &; Bromley (2004), predicted the cumulative size 



distributions for sub-kilometer objects to be within the ranges of the three models mentioned above 
and hence are not shown here. 

We next calculate the number of events for three models assuming 300 clear nights per year, 
eight hours of observation per night and 60 guide stars with 0.02 mas and SNR=100 for PSl three- 



year lifetime as show in Figure 11 From the TAOS upper limit, the most optimistic prediction for 
the number of events is ~100. The P&:S05 model with different break diameter have number of 
events ranging from a few to a few tens. For B04, we expect less than one event in PSl's three 
year life-time. 



7. Guide Stars Catalog 

As we have shown in the previous section, the detection efficiency is higher for smaller stellar 
angular size and higher SNR. Therefore, it is crucial to select guide stars with the smallest possible 
angular size yet with good SNR in order to maximize the detection rates. Nevertheless, it is 
not clear whether a star with SNR=120, 0^=0.07 mas is a better choice than another star with 
SNR=100 and 6^,=0.03 mas. To answer this question we chose a size distribution from models, such 
as P&;S05 with break diameter Dy, = 40 km, and calculated the number of expected events for a 



given star. In Figure 12 we show the total number of expected events per star in PSl's three- year 
lifetime as a function of angular size and SNR. Using this relation, we can now estimate the event 
rates for all matched stars (i.e., candidate guide stars) from Tycho2 and 2MASS catalogs and then 
select stars that give the highest event rates as our guide stars. 



In Figure 13 we show the distribution of my < 13 mag for all Tycho2 stars visible to PSl 3tt 
sky along with all matched stars between Tycho2 and 2MASS catalogs. For this work, we removed 
all matched stars with my < 8.5 mag due to possible saturation. As it can be seen, most of the stars 
with my < 12 mag were matched. On average, there are about 180 candidates with my < 11.5 
mag for every PSl 7 deg^ field. Furthermore, Figure 12 shows the number distribution of matched 



stars as a function of SNR and 0^. The distribution of SNR peaks around 80 and the angular size 
peaks around 0.02 mas. 

Considering that only 60 guide stars are needed for each target field, we selected 3 x 10^ guide 
stars from a total of ~ 1.27 X 10^ candidates!^ The dotted line in Figure! 131 is the distribution of 



®For 3 TT sky and PSl 7 deg^ field of view tliere are about 5000 unique target fields. 

''We have compiled a complete list for all matched stars with rriv < 13 mag. For the complete list see 
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(V-K) (V - K) from CADARS predicted angular size (mas) 



Fig. 9. — : Left panel: Distribution of {V — K) of all matched stars between Tycho2 and 2MASS 
catalogs. As it can be seen, there are two peaks in the distribution, indicating contribution from 
main sequence (left peak) and giant (right peak). Middle panel: The (V — K) distribution of the 
measurements from CADARS. The number of main sequence stars dropped quickly for (V — K) > 
2.0 and so is the number of giants for iy — K) < 2.0. Right panel: The distribution of predicted 
angular size of the matched stars between Tycho2 and 2MASS. Main sequence stars peak around 
0.02 mas, and giants peak around 0.05 mas. The overall distribution (solid-line) peaks around 
0.02 mas. 




Fig. 10. — : The straight solid line (P&S05) has q = fy. The slope after the break diameter has 
g = 3. Double power law from Bernstein has al = 0.88 and a2 = 0.32 (B04). TAOS (Z08) upper 
limit was anchored at 28 km on the B04 model and set an upper limit with q < 4.6 at 0.5 km. The 
arrow indicates an upper limit. 
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Fig. 11. — : Total number of expected events in the PSl three-year Hfetime, assuming all 60 guide 
stars are 0.02 mas and SNR=100, for different models. The arrow between Z08 and B04 indicates 
that the number of expected events will move along the dash-line as the upper limit (i.e. q) becomes 
smaller (as stronger constraints apply to the upper limit). Number of expected events from Z08 
gives the most optimistic value of 100. The number of expected events from P&S05's models 
range from a few to a few tens. For B04, we expect less than one event in the PSl three-year 
lifetime. 
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Fig. 12. — : Total number of expected events in the PSl three-year lifetime. Given SNR, 6^ and 
Pan & Sari's model with 40 km break diameter with q = 3, assuming all 60 guide stars have the 
same angular sizes and SNRs. 
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the top 3 X 10^ stars that were selected for the highest event rates. It is clear that most of these 
stars are also bright {my < 11.2 mag). In reality, even if the guide stars were selected based solely 
on their SNR, only a few stars per field would have not been selected from our selection based on 
the event rates. For each PSl target field, we provide multiple choices of guide stars ranked by the 
estimated event rates. Multiple choices will ensure freedom of selecting guide stars away from gaps 
in CCD arrays or dead pixels that would develop over time. 



8. Engineering Data 



A set of video mode engineering data (~ 166 stars-hours) was acquired in fall of 2008. The 
guide stars used in these images were randomly picked bright stars as opposed to being selected 



based on the optimized event rates as we described in previous section. Each FITS (Wells et al. 



1981 1 cube in the data set contains a few thousand 100 x 50 pixels images with one bright guide 



star. For the purpose of testing, these images were taken with 10 Hz instead of 30 Hz sampling. 
Nevertheless, we calculated all the related threshold values for 10 Hz sampling and tried to look 
for possible events in these engineering data. 



All images in the FITS cube were processed with SExtractor (Berlin Sz Arnouts 19961 to 



extract lightcurves for every guide star. Lightcurves for stars on the edges were excluded, and the 
first and last 100 data points in the lightcurves were removed because no stars were in these images 
due to shutter operations. 

We applied high-pass filter to the lightcurves to remove the low frequency components (cal- 
culated by a running mean from 21 data points) and then de-trending to remove the common 



features existing among the lightcurves taken at the same time (Kim et al. 20091. Figure 15 shows 



two examples of lightcurves along with their autocorrelations before and after the lightcurves were 
high-pass filtered and de-trended. The raw lightcurves showed considerable fluctuations due to ei- 
ther changing sky transparency or image motion, resulting in significantly lower measured SNR than 
predicted SNR from Equation[T] Measured SNR here is the median value divided by the standard 
deviation of the lightcurve. Also, strong autocorrelations were present in the raw lightcurves, which 
can compromise the detection technique. Nevertheless, the filtered and de-trended lightcurves show 
improved SNR and low autocorrelations. Dashed lines in the right panels are the two-sigma range 
for autocorrelation of white noise: (z — l)/z (z + 2), where z is the number of data points in the 



lightcurve and / is the lag in the autocorrelation (Box Sz Pierce 19701 



Finally we searched for possible events using the detection algorithm described in Q To 
determine if the performance of detection method on the engineering data is consistent with the 
simulations. One of the guide star (a = 333.4820917°, 6 = 1.2981389°) hghtcurves with ~ 
0.04 mas and my = 11.12 mag and measured SNR = 76.3 was selected and events with various 



http: / /timemachine. iic.harvard.edu/andrew/ps.matched.all.uniq 
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Fig. 13. — : The left panel: Histograms of all stars visible to PSl Stt sky in Tycho2 catalog with 
mv < 13 mag (TY2), cross-matched stars between Tycho2 and 2MASS (TY2-TMC) and the top 
300 thousands cross-matched stars (300K) that give the highest event rates based on their predicted 
angular sizes, SNRs and P&SOS model. The peak distribution is near my ~ 10.8 mag which has 
predicted SNR ~ 100. The right panel: Cumulative distribution. On average, there are 180 stars 
with mv < 11.5 mag for each PSl 7 deg^ field; hence we can have multiple choices of guide stars 
ranked by their predicted event rates. 
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sizes were injected. The detection method was then used to recover the injected events and therefore 
the efficiency (see Figure[T4]). Compared with the efficiency of 30 Hz for the similar SNR and 0*, 
one can see the detection efficiency of 30 Hz is better than that of 10 Hz. This also confirms that 
higher sampling rate has better detection efficiency. Meanwhile, the efficiency from the simulation 
for similar SNR and 0* was slightly better than the video mode lightcurve as expected. 

To further evaluate the quality of the lightcurves, we compared the predicted SNR with the 
SNR of lightcurves before and after filtering and de-trending (see left panel of Figure[T6]) . As it can 
be seen in the plot, the measured SNR is never as good as the predicted ones. In fact, the SNR 
of the lightcurves of bright guide stars never exceed ~ 200, which suggests that the SNR is limited 



by scintillation noise (Young 1967 Gilliland Sz Brown 1992 1. The scintillation noise is estimated 



by Equation 11 in GilHland & Brown (1992): 



SNR- 



0.09^-2/3(^)1-75 exp(-/i//io)(2At) 



-1/2 



(10) 



where A is the diameter of the telescope in centimeters, x is the air mass, h = 3000 m is the 
altitude at Haleakala, ho = 8000 m is the scale height and At is the exposure time in seconds. 
The scintillation limit based on PSl parameters are shown as dotted line in the right panel of 



Figure 16 The engineering data has SNR < 200 which is consistent with the scintillation noise 
with PSl parameters. The immediate impact of the reduction in the measured SNRs is the number 
of expected events. 



In Figure 11 , we estimated the event rates from the distribution peak of predicted SNR(=100) 
and 6-1, = 0.02 mas of the guide stars candidates. If the distribution peak of measured SNRs is 
reduced from 100 to 50, then the number of expected events will be halved (estimated from SNR = 
50, 9i, = 0.02 mas in Figure 12). Nevertheless, given the wide range of expected events depending 
on the models and other assumptions, we are still able to impose constraints and eliminate models. 

We realized that these engineering data might not be taken from the best sky condition with 
optimal optical and mechanical performances. Significant discrepancy can be seen between the 



measured SNR and scintillation noise limit (Figure 16). Possible sources of the noise are image 
motion, defocusing and other unknown electronic noise. We expect the quality of the video mode 
images will be improved in the future, and better SNRs can be achieved as of optical and mechanical 
problems are solved. In the end, no event in the engineering data set passed the one false positive 
threshold as was expected from our event rate calculations. 
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7) (km) 

Fig. 14. — : The detection efficiency of a guide star in the 10 Hz engineering data (soHd square) 
with 9-1, ^ 0.04 mas, SNR ~ 76.3. Events with various sizes were injected in the Hghtcurve to 

check the detection efficiency. Compared the detection efficiency of 30 Hz, it is clear that higher 
sampling rate has better detection efficiency. Meanwhile, the efficiency from the 10 Hz simulation 
(solid circle) for similar SNR and 9^, was slightly better than the video mode lightcurve as expected. 
Most importantly, it also confirms that our detection algorithm works for PSl lightcurve and can 
be used to search for TNO occultations. 
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Fig. 15. — : These two lightcurves and their autocorrelations are from the same video mode batch. 
Left panel: The crosses are the raw lightcurves and the solid lines are the high-pass filtered and 
de-trended ones. High-pass filtering was done by removing a running mean of neighboring 21 data 
points from the lightcurves. Right panel: The autocorrelations of the lightcurves before and after 



processing. The dashed lines are the two-sigma range for autocorrelation of white noise (Box &; 



Pierce 19701. As it can be seen, high-pass filtering and de-tredning improved the SNR of the 



lightcurves and also removed the correlation from the lightcurves. 
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Fig. 16. — : Left panel: The predicted SNR versus the measured SNR before and after filtering and 
de-trending. The dotted line has slope = 1 to guide the eye. Even after filtering and de-trending, 
most of the measured SNRs were worse than the predicted ones. Right panel: The air masses 
versus the SNR~^. G&B92 is the scintillation noise limit described by Gilliland Sz Brown (1992|. 
The solid line is a fitting with fixed slope = 1.75 of all measured SNR. The significant discrepancy 
suggests that other sources of noise, such as image motion, defocusing and unknown electronic 
noise, also contributed to the measured SNR. 
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9. Result and Discussion 

We have made a pre-survey study of using lightcurves from guide star video mode images to 
search for occultation by TNOs near 43 AU. Simulations were made to calculate the null hypothesis 
distribution (Figure[6]) and detection efficiency for various angular sizes and SNRs. Under 30 Hz 
sampling, PSl can detect objects as small as ~ 400 m at 43 AU (Figure[7]). Using (V — K) 
indices, we predicted the angular sizes of all matched stars between Tycho2 and 2MASS catalogs 
with my < 13 mag and mx < 16 mag above the PSl southern declination limit 5 > —30°. The 
distribution of stellar angular sizes peak around 0.02 mas in PSl Sir sky (Figure[9]). On average, 
there are about 180 stars with my < 11.5 mag for PSl 7 deg^ field of view. Given that we can 
choose only 60 guide stars for each PSl target field, we ranked our guide star candidates by their 
number of expected events based on their angular sizes, SNRs and model by P&:S05. As mentioned 
before, multiple choices of guide star will ensure us the freedom of selecting guide stars away from 
gaps in the CCD array or dead pixels that would develop over time. Based on the differential size 
distribution from different models and the threshold to have one false positive in the PSl three- 
year lifetime, we estimated the number of expected events could be somewhere from 1 to ~ 100 



(Figure 11 1 . The engineering data allowed us to investigate the quality of the lightcurve and develop 
the detection pipeline for the upcoming real data. We have established that the detection technique 
performs as well with the filtered engineering data. We also realize that the true SNR is limited 
by systematics. Some of the systematics will be removed as we move to real data stream, however 
some of the systematics will remain. For example, scintillation will limit the SNR to about 200. We 
have recalculated the event rates with the worst case scenario and found that the event rates were 
compromised by a factor of four. Even with this pessimistic estimation the event rate that PSl will 
find can allow us to place a constraint on the size distribution and hence the evolution of the TNOs. 
Using the available engineering data and detection efficiency at 0.4 and 0.5 km, we were able to 
derive an effective solid angle ~ 2.3 x 10~^^ and ~ 1.2 x lO"^'' deg^ and set the 95% confidence 
upper limit on surface number density at A^(> 0.4 km) 1.3 x 10^^ and N{> 0.5 km) ~ 2.47 x 10^° 
deg-2 (Fi gure|17p. In future work, for the upcoming real data, we will set the threshold based 
on maximum true positive to false positive ratio, which allows more candidate events for further 
investigations. Meanwhile, the video mode lightcurves can also be used to search for objects in 
Sedna like orbits from 100 to 1000 AU. We will work on a new detection algorithm that is capable 
of searching for objects in this region. 



Work at NCU was supported by the grant NSC 96-21 12-M-008-024-MY3. Work at the CfA was 
supported in part by the NSF under grant AST-0501681 and by NASA under grant NNG04G113G. 
Work at ASIAA was supported in part by the thematic research program AS-88-TP-A02. And we 
also like to thanks IIC for hosting J-H Wang for his staying. 



- 25 - 



REFERENCES 

Allen, R. L., Bernstein, G. M., & Malhotra, R. 2001, ApJL, 549, L241 
Bailey, M. E. 1976, Nature, 259, 290 
Bernstein, G. M., et al. 2004, AJ, 128, 1364 
Bertin, E., & Arnouts, S. 1996, 117, 393 
Bianco, F. B., et al. 2009, ArXiv e-prints 

Bickerton, S. J., Kavelaars, J. J., & Welch, D. L. 2008, AJ, 135, 1039 

Bickerton, S. J., Welch, D. L., & Kavelaars, J. J. 2009, AJ, 137, 4270 

Box, G., &L Pierce, D. 1970, Journal of the American Statistical Association, 1509 

Chiang, E. I., & Brown, M. E. 1999, AJ, 118, 1411 

Fuentes, C. I., & Holman, M. J. 2008, AJ, 136, 83 

Gilliland, R. L., & Brown, T. M. 1992, PASP, 104, 582 

Gladman, B., et al. 2002, Icarus, 157, 269 

Gumbel, E. J. 1943, Ann. Math. Statistics, 14, 253 

H0g, E., et al. 2000, A&A, 355, L27 

Jewitt, D., & Luu, J. 1993, Nature, 362, 730 

Kaiser, N., & Pan-STARRS Team. 2002, in BAAS, Vol. 34, 1304 

Kenyon, S. J., & Bromley, B. C. 2004, Nature, 432, 598 

Kim, D.-W., et al. 2009, MNRAS, 770 

Lehner, M. J., et al. 2009, PASP, 121, 138 

Luu, J., et al. 1997, Nature, 387, 573 

Maiz-Apellaniz, J., et al. 2004, ApJS, 151, 103 

Nihei, T. C, et al. 2007, AJ, 134, 1596 

Pan, M., & Sari, R. 2005, Icarus, 173, 342 

Pasinetti Fracassini, L. E., et al. 2001, A&A, 367, 521 

Reed, B. C. 2003, Aj, 125, 2531 



- 26 - 



Rice, J. 2007, Mathmatical Statistics and Data Analysis, 97 

Roques, F., & Moncuquet, M. 2000, Icarus, 147, 530 

Roques, F., Moncuquet, M., k Sicardy, B. 1987, AJ, 93, 1549 

Roques, F., et al. 2006, AJ, 132, 819 

Russell, H. N. 1916, ApJ, 43, 173 

Scholz, M. 1997, in lAU Symposium, Vol. 189, 51-58 

Skrutskie, M. F., et al. 2006, AJ, 131, 1163 

Tonry, J., Burke, B. E., & Schechter, R L. 1997, PASP, 109, 1154 
Trujillo, C. A., et al. 2001, AJ, 122, 2740 
Van Belle, G. T. 1999, PASP, 111, 1515 

Wells, D. C, Greisen, E. W., & Harten, R. H. 1981, A&AS, 44, 363 

Young, A. T. 1967, AJ, 72, 747 

Zhang, Z.-W., et al. 2008, ApJL, 685, L157 



-27- 



B04 
Z08 
BKW 



0.4 



0.5 0.7 
D (km) 



1.0 



Fig. 17. — : The upper limit sets by PSl 22 star-hours engineering data at 95% confidence limit 
along with upper hmits by TAOS, BKW ([Bickerton et al.||2008|) and B09 ([Bianco et a"L]|2009[). 



